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Abstract 

An approximate solution technique has been de- 
veloped for three-dimensional, inviscid, hypersonic 
flows. The method uses Maslen’s explicit pressure 
equation and the assumption of approximate stream 
surfaces in the shock layer. This approximation 
represents a simplification of Maslen’s asymmetric 
method. The present method presents a tractable 
procedure for computing the inviscid flow over three- 
dimensional surfaces at angle of attack. The so- 
lution procedure involves iteratively changing the 
shock shape in the subsonic-transonic region until the 
correct body shape is obtained. Beyond this region, 
the shock surface is determined by using a marching 
procedure. 

Results are presented herein for a spherically 
blunted cone, a paraboloid, and an elliptical cone 
at angle of attack. The calculated surface pressures 
are compared with experimental data and finite- 
difference solutions of the Euler equations. Shock 
shapes and profiles of pressure are also examined. 
Comparisons of the results of the present method 
with experimental data and detailed predictions are 
very good. Since the present method provides a very 
rapid computational procedure, it can be used for 
parametric or preliminary design applications. A 
useful application would be to incorporate a heating 
procedure for aerothermal studies. 

Introduction 

An approximate solution of the three-dimensional 
(3-D) inviscid flow-field equations for hypersonic 
speeds is an important design tool, especially for 
the prediction of surface heating rates on reentry 
vehicles. Coupled inviscid- viscous calculations have 
been demonstrated to adequately predict the heating 
over a wide range of aerothermal environments. (See 
refs. 1 to 3.) The 3-D inviscid solutions used by these 
engineering aeroheating methods differ in complex- 
ity from simple modified Newtonian theory (ref. 1) 
to more exact finite-difference solutions of the 3-D 
Euler equations (refs. 3 and 4). Modified Newtonian 
theory is simple to use, but its accuracy is limited. 
Finite-difference solutions of the inviscid equations 
(refs. 4 to 7) may be too computer-intensive for a 
preliminary design study. A rapid but accurate 3-D 
inviscid analysis would significantly enhance current 
engineering aerothermal capabilities. 

Surface heating rates on reentry vehicles may 
also be obtained by numerically solving the full 3-D 
Navier-Stokes equations. (See ref. 8.) Solutions to 
the various subsets of the Navier-Stokes equations, 
including the parabolized Navier-Stokes (ref. 9) and 
viscous shock-layer (refs. 10 to 12) equations, have 


also been obtained for 3-D flows. Although these 
methods yield good results, they require excessive 
computer storage and CPU (central processing unit) 
time to be practical for a preliminary-design envi- 
ronment in which a range of geometries and flow 
parameters are to be studied. An approximate 3-D 
inviscid method coupled with an engineering aero- 
heating prediction method can easily be run on a 
computer workstation and complements the more ac- 
curate CFD (computational fluid dynamics) methods 
in the design of reentry vehicles. 

Maslen (ref. 13) developed an approximate invis- 
cid flow-field method for axisymmetric bodies in the 
early 1960’s. From simplifications made for hyper- 
sonic speeds, Maslen obtained an explicit expression 
for the first-order pressure along lines normal to the 
shock wave. Since the body shape is determined from 
a given shock shape, this method is referred to as an 
inverse technique. Solving for a specified body shape 
can be accomplished indirectly by iteration of the 
shock shape. (See refs. 14 and 15.) This technique 
gives reasonable results, but the first-order pressure 
equation is inaccurate on the conical afterbody of a 
spherically blunted cone. 

Maslen subsequently outlined a method for three- 
dimensional inviscid flow (ref. 16), and improved the 
explicit pressure equation by deriving a second-order 
approximation of the momentum equation normal 
to the shock. Although more accurate, the 3-D 
method is complicated and unwieldy. In fact, only 
preliminary results were presented for blunt-body 
flows. (See refs. 17 and 18.) 

The present analysis simplifies the asymmetric 
problem by using two stream functions that approx- 
imate the actual stream surfaces in the shock layer 
and a reduced form of the Maslen second-order pres- 
sure equation. Thus, the present method, while ap- 
plicable to three-dimensional flows, retains the util- 
ity and accuracy of Maslen’s axisymmetric technique. 
Comparisons of the present results in terms of pres- 
sure profiles, surface distributions, and correspond- 
ing shock shapes are presented with Euler predictions 
and experimental data for spherically blunted cones, 
ellipsoids, and paraboloids at angle of attack. 

Symbols 

f?, 6, c, d conic section parameters 

(eqs. (19) and (20)) 

e^, e r , unit vectors of cylindrical 

coordinate system 

e^, e^, e n unit vectors of curvilinear 
coordinate system 



/ 

shock radius (eq. (1)), non- 
dimensionalized by L 

9 

metrics of curvilinear coordi- 
nate system (see appendix A) 

> hp > bn 

scale factors of curvilinear co- 
ordinate system, nondimen- 
sionalized by L 

J 

Jacobian of coordinate 
transformation 

L 

reference length 

M 

Mach number 

n 

coordinate in normal direction, 
nondimensionalized by L 

P 

position vector from origin of 
coordinate system 

V 

static pressure, nondimension- 
alized by pooV 2 ) 

R 

radius of curvature, nondimen- 
sionalized by L 

r 

radius in cylindrical coordinate 
system, nondimensionalized 
by L 

u, v, w 

velocity components of curvi- 
linear coordinate system, 
nondimensionalized by 

V 

velocity magnitude 

V 

velocity vector, nondimension- 
alized by Vqo 

x,y,z 

Cartesian coordinate system 
referenced to shock or body, 
nondimensionalized by L 

a 

angle of attack 

M 

coordinate perpendicular to 
shockline and unit normal 
vector, nondimensionalized 
by L 

r 

shock-wave angle relative to 
free-stream velocity (see fig. 1) 


shock- wave angle in circumfer- 
ential direction (see fig. 1) 

e 

axial distance from origin at 
which integration of shock 
variables is started 


rj ratio of stream function ip to 

its value on shock 

6 inclination angle of body 

surface with respect to body 
axis 

£,£ curvilinear coordinates along 

shock line, nondimensionalized 
by L 

p density, nondimensionalized by 

Poo 

a shock-wave angle, <p — 6$ 

$,'0 stream functions (eq. (9)), 

nondimensionalized by 

PooVooL 2 

p cylindrical coordinate (see 

i) 

Subscripts: 

b value on body 

c property on conical afterbody 

of spherically blunted cone 

i conic-section parameters 

specific to p = Constant plane 

j marching step 

s value on shock surface 

oo * free-stream conditions 

Analysis 

Coordinate Systems 

For a three-dimensional shock wave, a cylindrical 
coordinate system (x,r, <f>) is used, and the rr-axis 
is aligned with the free-stream velocity vector. The 
shock- wave geometry (ref. 19) may be described by 


r s = f(x, 4>) 

(i) 

Two angles that define the shock- wave 

shape are 

, c 1 df 

tmS *= m 

(2) 

n p 

tanT = ~ cos 8^ 

(3) 

A third angle is given by the relation 


, ~e- 
1 

-e- 

II 

b 

(4) 


All three angles are shown in figure 1. For the special 
case of axisymmetric flow, T = r(x), 6$ = 0, and 

(7 — <fi. 
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Next a shock-oriented curvilinear coordinate 
system (£,/?, n) is defined with corresponding unit 
vectors (e^,e^,e n ) and velocities (u,w,v). This co- 
ordinate system (refs. 1 and 20) is appropriate for 
hypersonic flow (M 0 0 » 1) and thin shock layers. 
The unit vector e n5 is the inward vector normal to 
the shock surface. The coordinates £ and (3 are cho- 
sen such that e^ s is in the direction of the tangen- 
tial velocity component just inside the shock surface. 
The unit vector ep $ is perpendicular to e^ s and e n5 . 
(See fig. 2.) This coordinate system is defined as or- 
thogonal at the shock surface but is nonorthogonal in 
the shock layer for a general three-dimensional shock. 
The derivations of the unit vectors and the metrics 
in the shock layer are outlined in appendix A. 

The postshock relations can be written as 


respectively. The scale factors h^ si hp s , and h ns gov- 
ern the stretching of the corresponding coordinates. 
The scale factor in the n-direction h ns is unity, since 
it is a straight-line distance. 

Governing Equations 

Stream functions . The continuity equation for 
three-dimensional flows given by 

V • (pV) = 0 (8) 

is satisfied by two stream functions. The continuity 
equation can be written as 

pV = Vi/> x (9) 


Vs = u s e^ s + v $ e ns (5) 

w s = 0 


whereas in general, through the shock layer, 

V = ue ^ + ve n + wep (6) 


The component of the free-stream velocity tangent to 
the shock wave is unchanged across the shock wave; 
therefore, the unit vectors e^, e^, and e n at the shock 
are 


e^ s = cos Te x + sin T (cos fye r 
ep s sin cos 


sin <5^e^) ^ 


} (?) 


e ns — sin Te*: - cos V (cos fye r - sin J 


where tp and $ represent the two stream functions 
(ref. 16). 

As noted in references 13 and 16 for blunt bodies 
at hypersonic speeds, most of the mass flow is near 
the shock wave, where the velocity component w 
is small. A simplifying assumption is made in this 
study that w = 0 throughout the shock layer. Thus, 
if <£ is set equal to /3, ip becomes 


d'lp _ puJ 
dn h £ 


( 10 ) 



( 11 ) 


The differential arc lengths along each coordinate 
direction at the shock are h^ s d£, hp s d(3 , and h ns dn, 


where J represents the Jacobian of transformation 
from (£, /?, n) to (x, y, z) and is given by 


J — h'£ ) sh'f3s 



r 1 9(7 \ 
ncos T — — - 

hfis dfi J 


V fyfc 9 ( 3 ) 


21 


( 12 ) 


These definitions of 'ip and are not unique, and they satisfy the exact flow-field equations only at the shock 
wave. The intersections of j3 = Constant planes with the shock surface are referred to as shock lines. Shock 
lines are in the direction of the ^-coordinate, and for axisymmetric flow, /3 = Constant planes are meridional 
planes. 

Pressure equation . The momentum equations for steady, inviscid flow may be written as 


p(V* VV) + Vp = 0 


(13) 
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Writing these equations in the £, /3, n system gives the following momentum equation normal to the shock: 


u dv dv 
h^d^ V dn 


1 - n^r -a- - n 


R R 


J_dF 

h/3s dP 


+ l^=o 

p on 


where R is the radius of curvature of the shock in the £-n plane. Transforming this equation from the £, /?, n 
system to a new set of independent variables (£,/?, 77 ) by using equations (10) and (11) gives 


dp _ / ip s I dv 
dp J | <9£ 


i-n-M i 


1 ar 

hp 8 df3 


h'Zsh/Ss't 7sinr 
J dp 


where 


/? = /? 
r\ - 

Note that 77 = 1 on the shock and rj = 0 on the body surface. Prom the transformation, an expression for the 
velocity component v normal to the shock can be obtained and is given by 

u dn , h is hp s ri sinT 

v= h,T ( + ^j — (16) 


The only assumption used in equations (15) and (16) is that the velocity component w is equal to zero 
throughout the shock layer. However, to obtain explicit expressions for the pressure and normal velocity 
component, additional assumptions are required. The following approximations, which are consistent with the 
simplifications in references 13 and 16, are valid for hypersonic flow (Moo >> 1) and thin shock layers: 


;) , ( ”- 1) 


P~ Ps 




Equations (15) and (16) can now be simplified as 


fcn\ (c , u s^s , i/’s tan T / 1 1 da\ r/ 2 - 1 
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and 


v(€,P,v) = v s (€,p) 


1 + 


0 s 

hp s cosT 


1 ■ _ 1 da \ .. 

— + cosT — — fa - 1) 

R hp s d(3 I 


(18) 


If hp s = r s , equation (17) becomes Maslen’s second- 
order pressure equation for axisymmetric flow. (See 
ref. 16.) The pressure and normal component of 
velocity can now be found explicitly along a line 
normal to the shock surface, since all variables in 
equations (17) and (18) are evaluated on the shock. 
The surface pressure is easily found, since the stream 
function 0 and, thus, rj are zero on the body. Equa- 
tions (17) and (18) are approximate. 

Other relations . The energy equation for 
steady, adiabatic, inviscid flow reduces to the sim- 
ple relation that the total enthalpy of the flow is 
constant. Also, for inviscid equilibrium or frozen 
flow, the (postshock) entropy is constant along a 
streamline. With these two relations, the density 
and streamwise component of velocity u can be com- 
puted. (See fig. 3.) The distance from the shock 
surface to the body along a line normal to the shock 
is calculated by integrating equation (10) and noting 
that the Jacobian and scale factors are functions of 
the distance n as outlined in appendix A. 

These relations for pressure, density, the two ve- 
locity components, and shock-layer thickness outline 
an inverse method of solution in which the shock 
shape (not the body shape) is known. Therefore, the 
shock shape must be changed until the correct body 
shape is produced. The resulting iteration procedure 
is handled differently in each region of the flow. 

Method of Solution 

In the stagnation region of a blunt body travel- 
ing at hypersonic speeds, the flow is subsonic. Be- 
cause of the elliptic behavior of the flow-field equa- 
tions in this region, a marching scheme is not well 
posed. Thus, the complete shock shape for the entire 
subsonic-transonic region must be determined iter- 
atively. A marching procedure is then used down- 
stream of the subsonic-transonic region, where the 
flow is completely supersonic. 

Subsonic-transonic region . In this investiga- 
tion, the blunt-nose region of the body in a body- 
oriented coordinate system is represented by a longi- 
tudinal conic section with an elliptical cross section 
as 

B bVb + zl = 2 R b x b - b b x\ (19) 

where R & is the nose radius of the body in the x^-z^ 
plane relative to the principal shock radius of cur- 


vature in the x-z plane, b & determines the longitudi- 
nal shape of the body, and B ^ governs the ellipticity 
of the body cross section. Van Dyke and Gordon 
(ref. 21) suggested that a conic-section body shape 
produces a shock surface that can also be described 
by a conic section. This assumption is used in the 
present analysis. Since a three-dimensional shock 
needs to be specified, longitudinal conic sections are 
blended in the circumferential direction by using one 
or more elliptical arcs to produce the shock surface in 
the subsonic region. One advantage to this approach 
is that parameters are easily added by blending more 
conic sections. 

The equation of the longitudinal (0 = Constant) 
conic sections is given by 

r 2 s + biXg - 2 CiX s + dix s r s = 0 (20) 


where b , , q . and di are parameters local to a plane 
where 0 is constant in the wind-oriented system. At 
each ^-location on the shock, the radii that are de- 
termined from equation (20) in each meridional plane 
are fitted with one or more elliptical segments to pro- 
duce the circumferential variation of the shock ra- 
dius. The number of parameters (fy, q, and di) to 
calculate then depends on the number of planes. In 
this report, only three planes (0 = 0, 0 = 7r/2, and 
0 = 7r) are blended by using one cross-sectional el- 
lipse to give a shock surface in the nose region. This 
method is adequate for axisymmetric shapes at angle 
of attack and for bodies with elliptical cross sections 
at moderate angles of attack. However, for elliptical 
bodies at larger angles of attack, more parameters 
and more elliptical segments may be needed to accu- 
rately describe the shock surface. 

In the iteration procedure, the shock shape is first 
assumed to be equal to the body shape, and the 
conic-section parameters for the shock are chosen 
accordingly. Next, the shock standoff distance is 
computed from the limiting form of the shock-layer 
equations on the stagnation line. (See ref. 1.) Shock 
lines are then traced from the stagnation region 
to the sonic surface by integrating the following 
differential equations that govern the shock shape 
along a shock line (/? = Constant): 


dr s 

dx 


tan T cos 6 $ 


( 21 ) 
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d<Ps 

tan T sin 6^ 

(22) 

dx 

r s 

dx 

= Ji0 s tanT 

(23) 

dh/3s 

dx 

da 

= tanr w 

(24) 

dsinT 

1 

(25) 

dx 

“ ~R 


where x is the independent axial coordinate. The 
shock variables are integrated along each shock 
line with a variable-step-size, third-order, predictor- 
corrector, ordinary differential-equation solver. An- 
gles and derivatives are determined by blending the 
longitudinal conic sections. The transformation op- 
erators used in the development of these equations 
are derived in appendix B. Since these equations are 
indeterminate at the stagnation line (x s = 0), the 
integration begins a small distance e from the origin 
(ref. 19). At x s = e, the coordinate is assumed to 
be equal to the curvilinear coordinate ■/?. 

The approximate shock-layer equations for pres- 
sure, density, velocity, and shock-layer thickness are 
used to calculate the position of the body at dis- 
crete locations. These points are placed primarily 
near the end of the subsonic-transonic region, where 
a good starting solution is crucial to the downstream 
marching procedure. The error between the geomet- 
ric and calculated distance from the shock to the 
body is used to change the parameters of the shock 
surface in a quasi-Newton nonlinear equations solver. 
If the number of parameters is less than the number 
of calculated body points, the shock is determined 
in a least-squares sense. However, improved results 
downstream are obtained if the points on the body 
surface are matched exactly. Because three planes 
are used in this report, a total of nine parameters 
are varied until nine points on the body have the 
correct shock-layer thickness. 

Supersonic region . Once past the transonic 
region, the flow is totally supersonic and a marching 
scheme is well posed. The shock surface and resulting 
shock lines from the transonic region form a starting 
solution for the marching procedure. 

The differential equations that govern the shock 
shape along a shock line are given by equations (21) 
to (25). However, in the supersonic region, 6$ and 
da/d/3 are computed from finite differences in the 
circumferential direction at each x-location. 


The general marching procedure is outlined here. 
Starting at xj, predicted values for the shock vari- 
ables (r s ,0 s ,^ 5 ,/i^ 5 , and T s ) are obtained at xj+\ 
for each shock line. The angle 6$ and the deriva- 
tive da/d/3 are computed from finite differences by 
using the predicted values of the shock radius. On 
each shock line, the local shock curvature 1/R is var- 
ied in the approximate shock-layer equations that 
govern pressure, density, velocity, and shock-layer 
thickness until the calculated thickness matches the 
geometric shock-layer thickness. Convergence usu- 
ally required two or three iterations with the secant 
method. These values of the shock curvature are now 
used to determine corrected values of the shock vari- 
ables at Xj+i for each shock line. Updated values for 
8^, da/d/3, and 1/R are calculated in the same fash- 
ion as in the predictor step. Upon completion of the 
corrector step at Xj+\ for each shock line, predicted 
values for the shock variables are obtained at xj+ 2, 
and the entire process is repeated. 

Results and Discussion 

Results at perfect-gas conditions are presented for 
this investigation over a spherically blunted cone, a 
paraboloid, and an elliptical cone at angle of attack. 
Solutions are described in a body-oriented coordinate 
system with the circumferential angle equal to 0° 
in the windward plane of symmetry and 180° along 
the leeward ray. Distances are referenced to the nose 
radius of the body. Nondimensional pressures are 
obtained by dividing by PooV^. 

Surface pressures obtained with the present tech- 
nique for a spherically blunted circular cone at 0° 
angle of attack are shown in figures 4 and 5. The 
cone half-angle is 10° for both cases, and the free- 
stream Mach numbers are 5 and 20, respectively. 
Comparisons with a finite-difference solution (ref. 7) 
of the inviscid-flow equations show excellent agree- 
ment with the higher Mach number results of the 
present method. Comparisons of the detailed and 
present results at the lower Mach number condition 
are good. The discrepancies at this condition may be 
the results of shock-layer approximations used in the 
present method. 

The pressure distribution over the 10° blunted 
cone at an angle of attack of 10° is shown in fig- 
ure 6. The solution in the wind axes over the spher- 
ical nose at angle of attack is identical to the nose 
solution over a sphere at 0° angle of attack if the 
subsonic-transonic region remains on the spherical 
cap. The surface pressures along the windward ray 
again compare favorably with the predicted values 
of reference 7. However, in the leeward plane, the 
present technique underpredicts the Naval Surface 
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Weapons Center (NSWC) surface pressure by an av- 
erage of 10 percent. The circumferential pressure dis- 
tribution in figure 7, the shock shape in the planes 
of symmetry in figure 8, and the pressure profiles in 
figures 9 and 10 confirm the high degree of accuracy 
of the technique in the region of the windward plane. 

Shock shapes and surface pressures for a parab- 
oloid at 8° angle of attack are presented in figures 11 
to 13. Experimental results are presented at free- 
stream Mach numbers of 9.9 for the surface pressures 
and 5.73 for the shock shape. Although a paraboloid 
is axisymmetric, the shock shape produced in the 
nose region is fully three-dimensional when the body 
is at angle of attack. This comparison is a more sig- 
nificant test of the capabilities of the iteration pro- 
cedure in the subsonic-transonic region than for the 
spherically blunted cones at angle of attack, because 
in the wind axes, the flow is no longer axisymmetric. 
For the paraboloid at angle of attack, three longitu- 
dinal conic sections were blended by using a cross- 
sectional ellipse to produce the shock shape in the 
nose region. Six iterations using the quasi-Newton 
nonlinear equations solver were required for conver- 
gence. Good agreement (within 8 percent) in surface 
pressures between the present method and the exper- 
imental data (ref. 22) is shown in figures 11 and 12. 
In figure 13, the calculated shock lies slightly closer 
to the body than does the experimentally determined 
shock shape, but the agreement is good. Compar- 
isons in reference 22 at 0° angle of attack were made 
with the axisymmetric Maslen technique of Zoby and 
Graves (ref. 14), and a similar result was observed. 

Surface pressures over an elliptical cone with a 
cone half-angle of 10.26° in the windward plane and 
an ellipticity of 1.5 are shown in figures 14 and 15 
at a free-stream Mach number of 10. The experi- 
mental data (ref. 23) displayed are for a pointed el- 
liptical cone. Since the present method currently al- 
lows for blunted bodies only, a very small bluntness 
is assumed for the elliptical cone. Circumferential 
pressures are shown for a position far downstream 
(x/Rb = 120), where the surface pressures should ap- 
proach sharp cone values. As in the case of a spher- 
ically blunted cone at angle of attack, the surface 
pressures on the elliptical cone at 10° angle of attack 
show excellent agreement with the experimental data 
near the windward plane, but they diverge as the lee- 
ward ray is approached. However, viscous effects are 
more pronounced in the leeward plane, and an invis- 
cid method is not appropriate for calculations in this 
region. 

To estimate the relative computing times of the 
present technique and more exact methods that nu- 
merically solve the full inviscid equations, the present 
method and the Euler equation solver given in refer- 


ence 7 were run on a CDC® Cyber 170 computer and 
run times were compared. Solutions over spherically 
blunted cones with cone half-angles of 5° and 10° and 
lengths from 50 to 100 nose radii were computed for 
angles of attack of 0° to 10°. The Mach number for 
all cases was 10. Based on these comparisons, the 
present method is consistently five times faster per 
meridional plane than the Euler code of reference 7. 
For a 10° cone at 10° angle of attack, and for a length 
of 50 nose radii, the Euler code requires 150 CPU sec 
of computer time, while the present technique re- 
quires 30 sec. Nine meridional planes were used in 
the Euler code, and a corresponding number of shock 
lines were used in the present method. 

Solutions were also computed over the nose region 
of a blunt elliptical cone with an ellipticity of 1.5 us- 
ing the present method and the time-dependent Eu- 
ler solver code HALIS (High Alpha Inviscid Solution) 
of reference 4. Both methods were run on a Cray-2 S 
supercomputer, and run times were compared. For 
a length of 0.75 nose radii and a free-stream Mach 
number of 10, HALIS requires 120 CPU sec, while the 
present technique requires less than 1 CPU sec. Sur- 
face pressure distributions in the planes of symmetry 
are presented in figure 16. The circumferential pres- 
sure distribution at an axial location of 0.4 nose radii 
is shown in figure 17, and there is excellent agreement 
between the present method and HALIS. These com- 
parisons indicate that the present technique may be 
much faster than more exact CFD methods for 3-D 
nose shapes. 

The present approximate 3-D inviscid method has 
been shown to produce pressure profiles and surface 
distributions and corresponding shock shapes that 
are in good agreement with detailed predictions and 
experimental data. The procedure significantly im- 
proves the current engineering capability of predict- 
ing inviscid flow fields over 3-D surfaces. The present 
inviscid technique may be coupled with the aero- 
heating prediction method described in reference 1 
to calculate heating rates near the windward plane 
of symmetry on reentry vehicles. The speed of the 
present method allows it to be easily run on a com- 
puter workstation; therefore, it is ideal in a design 
environment. 

Concluding Remarks 

An approximate technique has been developed 
that determines the flow field over blunt- nose bodies 
at hypersonic speeds. An explicit expression from 
Maslen for the pressure across the shock layer and 
the assumption of approximate stream surfaces in the 
shock layer are used to simplify the three-dimensional 
flow. The method is applied to the solution over 
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spherically blunted cones, paraboloids, and elliptical 
cones at angle of attack for a perfect gas. 

The present technique predicts shock surfaces and 
shock-layer properties that compare favorably with 
the more exact numerical solutions and experimen- 
tal results. Excellent agreement is obtained near 
the windward plane of symmetry. The method is 
also very rapid when compared with the more exact 


finite-difference solutions of the inviscid equations, in 
particular the solution over three-dimensional (3-D) 
bodies. Therefore, this technique is especially attrac- 
tive for engineering aeroheating methods that require 
inviscid hypersonic flow-field information. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
July 11, 1990 


8 



Appendix A 


Variation of Unit Vectors and Metrics Across Shock Layer 


If P s is the position vector, relative to the origin, for a point on the shock wave and P is a position vector 
for a point normal to the shock wave at P s , then 


P(£, P, n) = P s (£, P) + / 3 ) 


(Al) 


The coordinate directions in the shock layer are found from the following equations: 



d P 

e s = 

d£ 


d P 

e /3 = 

dp 


d P 

— 

dn 


(A2) 


Differentiating equation (Al) and substituting into equation (A2), using the definition of the unit vectors at 
the shock (eq. (7)), yields 


e £ — 




, i\ i ar 

i _ n _j e ^ + n __ e/3s 

i dr 

hfis dP 


1 dr ( _ 1 do \ 

»^^e s ,+ (l-»cosr — 


(A3) 


Gn — Gn 


where 


± i_dr 

R ~ h^ s d£ 

1 dT _ cosT da 
h/3 s dP h$ s d£ 


The metrics for a nonorthogonal coordinate system are given by 


Qij — * G j 


(A4) 
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where % and j are indices that represent the coordinates f , /3, and n. Therefore, 


5« = 

5/3/3 = 
9nn — 1 


( i r 

90 - 9 ^ - hshps I n h ^ df3 


1 i \ 2 / i dr N 

1 — n— I + In- — — 

R) V h 0s dp. 


\ i da \ 2 ( 1 dr \ 2 

y ncos %,w + (%iw) 


(A5) 


1 - n i) + { 1 - naxT Kf t %, 


9^n — 9n £ — 0 
9pn ~ 9n/3 ~ 0 


The Jacobian of transformation from (£, /3,n) to (x,y,z) is given by 


J \9ij\ ^ 



1 da 


i or 


n cos T - — — - n- — — 


hp s dp 


h(3s dfl j 


(A6) 


The scale factors are defined to be the magnitudes of the vectors and are therefore 


h = (5«) 1/2 

h 0 = (5/3/3) 1/2 > 

hn ( 9nn — 1 , 


(A7) 
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Appendix B 

Transformation Operators at Shock Surface 

The unit vectors at the shock in the curvilinear coordinate system (£,/?, n) are 

e^ s = cosTe x + sin F (cos — sin^e^) ' 

<3{1s = sin ^e r + cos 6^ > 

e ns = sin T e x - cosr(cos^e r - sin _ 


(Bl) 


If the transformation operators are derived for the curvilinear system with respect to the cylindrical coordinates 
and a position vector P s on the shock surface is defined, the vector may be expressed as 


where 

The derivative of this vector is 


P.S = xe x + fe r 
r s = f{x, (p) 


(B2) 


dP s = ( e 


+ l e ') dl + ( e * + 7^ er ) /<w 


In curvilinear coordinates, equation (B3) becomes 

9P 3P 

dPs = d(, + ~ d/3 = e^ s + e [js h fis d(3 

Using equations (B3) and (B4) and taking the dot product with e x gives 

dx dx 

dx e^ s - G x h^ s + Gp s • G x hp s d(3 — d £ 4- dfi 


(B3) 


(B4) 


(B5) 


This process yields 


dx 

= e £s * e x^s 


dx 

Qp= e l3s- exhp s 

Taking the dot product of equations (B3) and (B4) with respect to <j> gives 

d(j> e £s * , 

9? = ~l~ h ’ 

d(f> e /3s ' 7 

9/r 

The transformation operators at the shock may now be written, using equations (Bl), (B6), and (B7) as 

Id _ d sin T sin 6a d 

— — = cos r — — 

h^ s dx r s d<p 

1 d cosS 6 d 


(B6) 


(B7) 


hp 8 d/3 r s dp 


(B8) 
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Shock 




Rear View 


Figure 1. Shock wave and body geometry. 




Figure 3. Coordinate positions (/3 = Constant) 


Present method 





Figure 5. Axial-surface pressure distribution on blunt cone with 6 C = 10°, a = 0°, and M ^ = 20. 



♦b-° 



180 ° 


Figure 6. Axial-surface pressure distribution on blunt cone with Q c = 10°, a = 10°, and M x = 10. 
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4> b . deg 

Figure 7. Circumferential surface pressure distribution on blunt cone with 6 C = 10°, a — 1( 
x b /R b = 9.6. 



x b /R b 


Figure 8. Shock shape in windward and leeward planes of symmetry for blunt cone with d ( 
and M™ = 10. 




Figure 10. Radial profile of pressure with 6 C = 10°, a = 10°, Moo = 10, x^/Rf, = 9.6, and </>j 









X b /R b 


Figure 13. 


Shock shape in windward and leeward planes of symmetry for paraboloid with a — 8° and 


Moo -5.73. 



Figure 14. Circumferential surface pressure distribution on 1.5 elliptical cone with a = 0°, M^ = 10, and 

e c - 10.26° {< t> b - 0°). 
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Figure 15. Circumferential surface pressure distribution on 1.5 elliptical cone with a = 10°, M 0 0 = 10, and 
9c = 10.26° {<t> b = 0°). 



Figure 16. Axial-surface pressure distribution on 1.5 elliptical nose with a — 0° and = 10. 
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